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ABSTRACT 


The  mechanics  of  quasi-static  crack  closing  and  bonding  of  surfaces  of 
the  same  or  different  linear  viscoelastic  materials  is  described.  Included  is 
a  study  of  time-dependent  joining  of  initially  curved  surfaces  under  the 
action  of  surface  forces  of  attraction  and  external  loading.  Emphasis  is  on 
the  use  of  continuum  mechanics  to  develop  equations  for  predicting  crack 
length  or  contact  size  as  a  function  of  time  for  relatively  general 
geometries;  atomic  and  molecular  processes  associated  with  the  healing  or 
bonding  process  are  taken  into  account  using  a  crack  tip  idealization  which  is 
similar  to  that  used  in  the  Barenblatt  method  for  fracture.  Starting  with  a 
previously  developed  correspondence  principle,  an  expression  is  derived  for 
the  rate  of  movement  of  the  edge  of  the  bonded  area.  The  effects  of  material 
time-dependence  and  the  stress  intensity  factor  are  quite  different  from  those 
for  crack  growth.  A  comparison  of  intrinsic  and  apparent  energies  of  fracture 
and  bonding  is  made,  and  criteria  are  given  for  determining  whether  or  not 
bonding  can  occur.  Examples  are  given  to  illustrate  use  of  the  basic  theory 
tor  predicting  healing  of  cracks  and  growth  of  contact  area  of  initially 
curved  surfaces.  Finally,  the  affect  of  bonding  time  on  joint  strength  is 
estimated  from  the  examples  on  contact  area  growth. 
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1.  Introduction 

In  1966  Williams  [1]  observed  that  "...  adhesive  failure  is  not  normally 
analyzed  as  a  fracture  problem  in  continuum  mechanics."  He  then  demonstrated 
that  cohesive  and  adhesive  fracture  may  be  analyzed  by  similar  methods  through 
considerat ion  of  the  stress  singularity  at  crack  tips  and  an  energy  balance. 
Following  this  pioneering  paper  the  understanding  of  adhesive  failure  has 
advanced  greatly,  clearly  aided  by  the  use  of  fracture  mechanics  theory  as 
urged  by  Wil liams. 

Quasi-static  analysis  of  adhesive  and  cohesive  fracture  of  linear 
viscoelastic  media  has  also  reached  a  relatively  mature  state  through 
application  of  fracture  mechanics  principles  le.g.,  2-7].  There  are,  however, 
still  important  unsolved  problems  of  cohesive  and  adhesive  fracture  of 
nonhomogeneous  media;  in  many  cases  the  stresses  on  the  crack  plane  are 
functions  or  the  local  viscoelastic  properties,  which  significantly 
complicates  the  analysis  [5,8].  Much  progress  has  been  made  as  well  on 
dynamic  fracture  of  linear  homogeneous  viscoelastic  materials  by  Walton  [9], 
among  others;  at  sufficiently  high  crack  speeds  the  stresses  depend  on 
viscoelastic  properties. 

The  problem  which  is  the  opposite  of  fracture,  the  bonding  of  surfaces, 
has  received  far  less  attention  in  the  context  of  fracture-like  mechanics. 
For  elastic  media  under  quasi-static  conditions,  the  analysis  of  bonding  is 
essentially  the  same  as  that  of  fracture.  Johnson  et  al .  [10]  and  Roberts  and 
Thomas  [11]  studied  by  experiment  and  linear  theory  the  bonding  and  unbonding 
of  flat  and  curved  surfaces  of  rubber  with  other  materials.  They  used  energy 
methods  analogous  to  those  employed  in  elastic  fracture  mechanics,  recognizing 
that,  contact  mechanics  and  fracture  mechanics  have  to  be  combined.  Although 
some  viscoelastic  effects  were  observed,  essentially  elasticity  theory  was 


2 


used.  Later,  Greenwood  and  Johnson  [12]  made  a  viscoelastic  fracture  analysis 
of  the  separation  of  a  sphere  from  a  plane,  allowing  for  dissimilar  materials, 
Put  assuming  there  are  no  shear  stresses  along  the  interface.  The  force  of 
attraction  acting  across  the  interface  was  based  on  a  theoretical  model  of  the 
force  between  planes  of  atoms  as  a  function  of  separation  distance. 

In  a  series  of  six  papers,  Anand  and  coworkers  [13-18]  studied  the  so- 
called  problem  of  autohesion,  which  is  the  bonding  of  surfaces  of  the  same 
material.  Contact  area  growth  between  irregular  (locally  curved)  surfaces  was 
analyzed  using  linear  viscoelastic  contact  mechanics.  The  influence  of  the 
force  of  attraction  between  surfaces  was  implicitly  neglected  compared  to  that 
due  to  a  constant  externally  applied  compression  when  predicting  contact  area 
growth.  The  increase  of  tensile  joint  strength  with  contact  time  was  assumed 
to  be  due  to  the  area  increase.  Other  studies  of  polymer-polymer  bonding  have 
assumed  complete  contact  is  achieved  immediately,  and  the  time-dependence  of 
strength  is  due  to  interdiffusion  of  molecular  segments  across  the  interface 
and  the  related  formation  of  entanglements;  many  references  may  be  found  in 
the  recent  paper  on  crack  healing  by  Kausch  et  al.  [19].  Jud  et  al .  [20]  and 
Stacer  and  Schreuder-Stacer  [21]  discussed  the  time-dependence  of  the  wetting 
or  contact  process  in  comparison  to  that  for  interdiffusion;  some  explicit 
estimates  of  contact  time-dependence  are  in  [20]. 

In  this  paper  we  address  the  problem  of  bonding  between  linear,  isotropic 
viscoelastic  media,  accounting  for  the  affect  of  interfacial  forces  of 
attraction  (or  surface  energy)  and  external  loading  on  the  rate  of  growth  of 
bonded  area.  Processes  are  assumed  to  occur  slowly  enough  that  inertia 
effects  are  negligible  and,  also  for  simplicity,  we  consider  only  the  mode  I 
problem  (in  which  only  normal  stresses  act  across  the  bond  surface).  The 
theory  shows  how  crack  healing  and  contact  area  growth  may  occur  under  either 
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external  tensile  or  compressive  forces,  and  a  relatively  simple  method  is 
developed  for  predicting  this  behavior. 

It  is  believed  the  results  will  be  useful  for  the  development  of  basic 
models  of  the  healing  of  damage  in  composites,  insofar  as  this  damage  consists 
of  micro-  or  macrocracks.  Significant  amounts  of  healing  have  been  observed 
experimentally  for  asphalt  [ 22 1  and  solid  propellant  [231,  which  are 
particulate  composites.  More  generally,  the  theory  may  be  used  to  account  for 
viscoelasticity  and  surface  force  or  energy  effects  in  adhesion  and  autohesion 
processes  when  externally  applied  compression  is  not  large  enough  to  produce 
complete,  immediate  surface-to-surface  contact.  Although  the  surface  energy 
may  be  small,  over  a  long  period  of  time  it  can  produce  considerable 
deformation  and  contact  area  growth  in  soft  materials.  With  a  minor 
modification  the  theory  in  this  paper  would  provide  a  means  for  predicting  the 
effect  of  a  concentrated  shear  stress  at  the  trailing  edge  of  so-called  "waves 
of  detachment"  observed  in  the  sliding  of  rubber  on  smooth  surfaces  124); 
addition  of  an  interfacial  shear  stress  may  be  done  in  the  same  way  as 
described  in  [5!  for  fracture. 

In  order  to  illustrate  the  type  of  behavior  studied  here,  consider  the 
problem  in  Figure  1.  Figure  la  shows  the  affect  of  (positive)  surface  energy 
on  the  contact  between  a  deformable  spherical  surface  and  a  rigid  flat 
surface.  Without  surface  energy.  Fig.  lb,  the  spherical  surface  is  tangent  to 
the  flat  surface  at  the  contact  edge,  and  the  load  necessary  to  produce 
contact  to  the  radius  ac  is  given  by  the  Hertz  theory  [25).  Less  compressive 
force  is  needed  in  Fig.  la  for  the  same  contact  radius;  indeed,  with  surface 
energy  a  tensile  force  can  be  supported.  The  dashed  lines  in  Fig.  1  indicate 
the  positions  of  undeformed  spherical  surfaces  after  experiencing  the  same 
amount  of  movement  at  the  remote  loading  point  as  in  the  deformed  body.  When 


a 


4 


the  body  is  viscoelastic,  the  contact  radius  varies  with  loading  rate  and,  at 
the  same  load,  may  be  considerably  different  for  positive  bonding  speed 
(ac>  0)  than  for  debonding  or  crack  growth  (a(.<  0)-  With  a  small  tensile  or 
zero  external  force,  the  interfacial  force  of  attraction  acting  across  the 
cusp-shaped  separation  at  the  contact  edge  causes  the  contact  radius  to  grow; 
if  the  long-time  modulus  of  the  sphere  is  zero  (e.g.,  noncrossl inked  rubber) 
the  contact  will  grow  indefinitely,  just  as  for  a  liquid.  The  force  required 
to  separate  totally  the  two  bodies  mey  be  large  if  considerable  viscoelastic 
dissipation  accompanies  the  separation  process.  The  separation  process  was 
analyzed  by  Greenwood  and  Johnson  (12]  using  the  author's  theory  of  crack 
growth  1 4 ] .  Here,  a  theory  is  developed  for  predicting  the  viscoelastic 
bonding  process  for  this  and  other  geometries  and  material  combinations,  as 
well  as  for  healing  of  internal  cracks. 

Section  2  introduces  the  notation  and  geometry  used  in  describing  the 
edge  of  the  bonded  area.  This  edge  region,  Fig.  2,  is  assumed  to  be  very 
small  compared  to  the  contact  size;  for  example,  it  exists  only  at  the  contact 
edge  in  Fig.  la  where  the  parabolic  surface  profile  touches  the  lower  flat 
plate,  without  the  local  tangential  separation  shown  in  Fig.  2,  the  stresses 
at  the  contact  edge  would  be  infinite,  according  to  linear  viscoelasticity 
theory.  With  edge  movement,  these  infinite  stresses  produce  infinite  strain 
rates  and  lead  to  physically  unacceptable  predictions  of  bonding  or  crack 
speed.  In  the  latter  case,  it  is  well-known  that  the  singular  stresses  lead 
to  a  physically  unacceptable  prediction  of  crack  growth.  Difficulties  arising 
from  use  of  the  classical  singular  stresses  were  recognized  by  Williams  |2], 
who  circumvented  the  problem  by  replacing  the  local  geometry  of  a  crack  by  a 
spherical  flaw.  Later,  Knauss  [3]  and  the  author  |4]  used  the  Barenblatt 
method  (or  what  is  often  called  the  Dugdale-Barenblatt  or  Dugdale  method)  to 


5 


remove  the  singularity.  The  latter  approach  is  used  in  this  paper;  it  amounts 
to  adjusting  8  and/or  ob  so  that  the  stresses  on  the  bond  surface  (and  in  the 
continua)  are  finite. 

Section  3  indicates  how  the  stresses  and  displacements  in  viscoelastic 
media  may  be  obtained  from  elastic  solutions.  This  approach  is  used  in 

Section  4  with  Barenblatt's  method  for  elastic  media  to  obtain  the  local 
viscoelastic  interfacial  stress  and  displacement.  In  contrast  to  the  crack 
growth  problem,  the  stresses  in  elastic  and  viscoelastic  bodies  are 

different.  This  difference  complicates  the  bonding  problem,  and  leads  to 
viscoelastic  effects  in  the  neighborhood  of  the  bond  edge  which  are 
considerably  different  than  predicted  with  crack  growth.  Nevertheless,  in 

certain  cases  it  is  possible  to  derive  relatively  simple  equations  for 
predicting  bonding  speed  and  the  bondinq-2one  length  8,  Fig.  2,  in  terms  of 
the  stress  intensity  factor.  In  Section  5  the  interface  bonding  stress 
ob.  Fig.  2,  is  assumed  constant,  and  exact  expressions  are  obtained  for  a 
material  characterized  by  a  generalized  power  law  creep  compliance; 

approximate  expressions  are  then  developed  for  a  more  general  form  of  creep 
compliance.  In  the  Appendix  a  simple  power  law  compliance  is  used  with  a 

bonding  stress  that  varies  with  distance  between  the  surfaces  to  derive  an 
exact  equation  for  bonding  speed;  the  dependence  of  bonding  speed  on  stress 
intensity  factor  is  found  to  be  the  same  as  for  constant  ob>  Section  6  is 
concerned  with  the  energy  needed  for  crack  growth  as  compared  to  that  for 

bonding  of  surfaces,  accounting  for  viscoelastic  dissipation  effects. 
Examples  of  crack  healing  and  the  bonding  of  initially  curved  surfaces  are 
given  in  Section  7,  using  the  theory  developed  in  Section  5.  The  theory  can 
be  readily  extended  to  certain  types  of  orthotropic  linear  and  nonlinear 
media,  as  discussed  in  the  concluding  remarks  in  Section  8. 
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2.  Description  of  the  contact  edge 

Referring  to  Fig.  2,  the  (x,y)  coordinate  system  is  stationary,  with  the 
y-axis  located  at  any  convenient  position  in  the  continuum;  y  =  0  defines  the 
surface  along  which  bonding  occurs.  We  assume  this  surface  is  locally  flat, 
but  do  not  restrict  its  shape  far  from  the  contact  edge  (relative  to 
b).  The  bodies  above  and  below  y  =  0  are  assumed  to  be  in  contact  to  the 
right  of  the  point  designated  as  P.  In  practice,  surface  irregularities  may 

prevent  complete  contact  of  the  adjacent  areas;  however,  for  our  purposes,  it 
will  be  sufficient  to  assume  only  that  the  scale  of  this  local  roughness  is 
small  compared  to  the  length  of  the  bonding  zone  e  which  is  defined  below. 

To  the  left  of  the  position  x  =  a-6  the  two  surfaces  are  assumed  to  be 
far  enough  apart  that  the  intermolecular  force  of  attraction  is  negligible. 
Over  the  length  a  the  force  of  attraction,  which  on  the  scale  of  the  continuum 
is  represented  by  the  tensile  bonding  force  per  unit  area,  ob,  tends  to  draw 
the  adjacent  surfaces  together  until  complete  contact  exists  (in  the  sense 
indicated  earl ier) . 

The  edge  of  complete  contact  (which  can  also  be  thought  of  as  a  crack 
tip),  whose  intersection  with  the  plane  of  the  page  is  the  point  P,  is  in 
general  a  curved  line  in  space  which  may  be  closed  or  else  open  and  ending  on 
the  edges  of  the  continuum.  It  is  necessary  to  assume  the  radius  of  curvature 
of  this  line  in  the  neighborhood  of  the  point  P  is  large  relative  to  6  in 
order  to  be  able  to  use  the  equations  of  plane  strain  in  a  local  bonding 
analysis;  this  neighborhood  is  defined  by  the  size  e. 

The  materials  above  and  below  the  bond  surface  are  assumed  to  be  linearly 
viscoelastic,  homogeneous  and  isotropic,  with  the  possible  exception  of  a  thin 
surface  layer  of  damaged  material;  this  layer,  which  may  be  the  result  of  the 
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surface  being  formed  in  the  first  place  by  means  of  crack  growth  or  any  other 
process,  will  have  no  practical  affect  on  our  analysis  if  its  thickness  is 
small  compared  to  a.  As  noted  previously,  the  scale  of  local  surface 
roughness  is  similarly  assumed  to  be  small.  These  surface  features  will,  of 
course,  affect  the  particular  values  of  ob  which  are  to  be  used  in  the 
theory.  Roughness  on  a  scale  much  larger  than  a  is  permitted. 

3.  Viscoelastic  stresses  and  displacements 

By  following  arguments  based  on  the  Laplace  transform  similar  to  those 
given  in  [4]  for  the  case  of  the  opening  mode  of  crack  growth,  we  can  express 
viscoelastic  stresses  and  displacements  in  terms  of  elastic  solutions  for 
crack  shortening;  this  approacn  amounts  to  the  use  of  the  extended 
correspondence  principle  for  contact  or  indentation  problems  (61. 

With  crack  growth  or  no  growth,  i.e.  a  >0  in  Fig.  ?,  the  normal  stress 
acting  along  the  plane  of  crack  prolongation  is  the  same  for  both  elastic  and 
viscoelastic  bodies  (which  are  assumed  to  be  identical  with  respect  to 
geometry  and  applied  loads)  if  this  stress  for  the  elastic  body  is  independent 
of  elastic  constants.  This  result  in  turn  leads  to  viscoelastic  displacements 
which  can  be  written  explicitly  in  terms  of  the  history  of  the  elastic  values. 

In  contrast,  it  may  be  shown  by  using  the  Laplace  transform  that  with 
crack  shortening  (or,  equivalently,  in  problems  with  growing  contact  area)  the 
normal  displacements  and  stresses  take  on  opposite  roles  with  respect  to 
elastic  solutions  if  the  elastic  normal  surface  displacement  is  independent  of 
the  elastic  constants.  This  latter  condition  is  not  necessarily  satisfied  for 
the  type  of  boundary  value  problems  of  interest  here  unless  the  Poisson's 
ratio  is  constant  and  all  traction  boundary  conditions  are  converted  to 
specified  displacement  boundary  conditions.  When  the  Poisson's  ratio  is 
constant  (or,  at  least,  the  effect  of  its  time-dependence  is  small  enough  that 
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it  may  be  assumed  constant)  the  more  straightforward  method  given  in  [261  for 
constructing  viscoelastic  solutions  from  elastic  solutions  may  be  used.  The 
method  is  not  based  on  Laplace  transform  theory  although  it  uses  a  so-called 
correspondence  principle,  designated  as  CP-Ill  in  (26).  This  correspondence 
principle  may  be  used  with  certain  types  of  nonlinearity  (including 
distributed  damage),  anisotropy,  and  aging.  However,  for  simplicity  we 
consider  here  only  linear,  isotropic,  nonaging  media. 

Let  ls  now  state  the  relevant  version  of  CP- III  and  then  define  the 
terms: 

If  sn./st  =  0  on  the  surface  Sy  and  dSy/dt  s  0  ,  the  viscoelastic 
solution  is 


•  =  'Ed; 

1J  1 J 


u  ; 


(1) 

(2) 


R  R 

where  and  u.  satisfy  the  field  equations  and  boundary  conditions  of  the 
so-called  reference  elastic  problem  (or,  more  briefly,  the  elastic  problem) 
defined  below. 

The  viscoelastic  and  elastic  bodies  are  assumed  to  be  geometrically 

identical  in  all  respects,  and  to  have  tractions  specified  over  the  portion  of 

the  (external  and  internal)  surfaces  designated  by  Sy  and  displacements 

specified  over  the  remaining  surfaces  Sy ;  the  complete  boundary  is  designated 

by  S.  Mixed  conditions  (e.g.  normal  displacement  and  shear  traction)  along  a 

boundary  may  be  specified;  but  for  rotational  simplicity  this  situation  is  not 

explicitly  stated  in  CP-III.  All  stresses  and  displacements  ui  in  the 

R  R 

viscoelastic  problem,  and  .  and  u.  in  the  elastic  problem  (where  i,j  = 
1,2,3)  are  referred  to  an  orthogonal  set  of  Cartesian  coordinates  x^  of  the 
undeformed  bodies.  The  n^  are  components  of  the  outer  unit  normal  at  any 
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N* 

V, 
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H  " 
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I 


point  on  the  undeformed  surface. 

The  braces  •;  }  denote  a  convolution  integral, 

f£df}  =  ER_1  A(t  -  t)  £  dx  (3) 

-  CD 

D 

where  f  is  a  function  of  time  such  as  or  o- j  ,  ER  is  an  arbitrarily 
selected  constant  (usually  with  units  of  modulus)  and  E(t)  is  the  uniaxial, 
linear  viscoelastic  relaxation  modulus;  for  an  elastic  material,  E  is  the 
Young's  modulus.  In  many  cases  the  quantity  (Ddf }  will  be  used;  it  is  defined 
using  the  creep  compliance  D(t), 


\ Ddf }  e  Er 


/oit-,) 


(4) 


«nere  D(t)  and  E(t)  are  interrelated  according  to 

;lE(t  -  O  §  d:  =  ;\(t  -  x)  37  dx  -  1  (5) 

o  o 

for  t  >  0;  also  E  =  0  =  0  for  t  <  0.  The  lower  limit  is  to  be  interepreted  as 
0’  in  order  to  include  the  singularity  in  the  derivatives  at  :=0. 

Additional  relationships  which  will  be  needed  later,  and  which  follow 
directly  from  (3)  -  (5),  are  concerned  with  inverses  of  the  linear  operators 
denoted  by  the  braces.  Namely,  if 


then 


g  =  -Ddf 


(6) 


f  =  <Edg;  (7)  i 

i 

j 

or,  equivalently, 


f  =  'Ed;0df;j 


and 


(8) 
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f  =  { Dd { Edf j )  (9) 

Equations  (8)  and  (9)  may  be  easily  verified  by  Laplace  transforming  them  and 
then  using  the  Laplace  transform  of  (5);  the  two-sided  Laplace  transform  can 
be  used  in  the  event  that  f  and  g  do  not  vanish  at'  negative  times. 

In  the  earlier  work  (261  we  assumed  that  all  stresses  and  strains  vanish 
for  t  <  0,  while  here  this  assumption  is  not  made.  One  can,  of  course,  always 
use  a  lower  limit  of  zero  in  (3)  and  (4)  by  shifting  the  time  axis  so  that  the 
body  is  unstressed  and  unstrained  for  t  <  0.  However,  when  using  CP- 1 1 1  in 
problems  for  which  dSy/dt  >  0  (such  as  crack  growth)  before  dSy/dt  <  0  (such 
as  crack  healing)  it  may  be  helpful  to  select  t=0  to  follow  the  dSy/dt  >  0 
process;  with  this  convention  stresses  and  strains  may  exist  at  t=0  without 
invalidating  CP- 1 1 1  if  they  can  be  calculated  without  regard  to  the  dSy/dt  >  0 
process.  Such  as  situation  arises  if  a  crack  first  grows  and  then  is  arrested 
for  a  long  enough  time  that  the  mechanical  state  of  the  continuum  at  t=0  is 
the  same  as  if  the  crack's  geometry  had  not  changed  when  t<0.  Consequently, 
one  can  solve  the  crack  shortening  problem  for  t  >  0  by  assuming  there  is  no 
change  in  crack  dimensions  for  t  <  0. 

Returning  to  CP-III,  the  elastic  problem  is  defined  by  the  standard 

P 

equilibrium  equations  with  a  body  force  F . , 

/3xj  *  fi  *  °-  <io> 

strain-displacement  equations, 

cjj  =  j(au*  /ax.  +  3U^  /3X.)  (11) 

and  stress-strain  equations, 


'y .'.y.v , 


11 


‘u  ■  '(1  *  »)0?J  -  -*,j  °kk  l/E„  (12) 

where  k  =  1,2,3.  The  usual  summation  convention  is  used  in  that  summation 
over  the  range  of  the  repeated  index  is  implied.  Also,  Er  and  v  are  the 
Young's  modulus  and  Poisson's  ratio,  respectively,  for  the  elastic  problem, 
and  6 .j j  is  the  Kronecker  delta  (i.e.,  <$.j  =  1  if  i  =  j  and  6 . j  =  0 
if  i  *  j).  The  boundary  conditions  are 


ij  J  i 


on  $. 


u-  =  U?  on  Su  (14) 

R  R 

The  quantities  F.  and  T.  specified  in  the  elastic  problem  are  not  the 
same  as  the  traction  T^  and  body  force  Fi  specified  in  the  viscoelastic 
problem  of  interest.  Rather,  they  are  to  be  found  using  the  convolution 
integral s, 


n  -  ;DdFi 


t?  s  ,-°dTi 


However,  the  specified  displacements  are  the  same  for  both  problems, 


U-  -  Ui  (17) 

The  viscoelastic  boundary  value  problem  is  defined  by  the  field  equations 
3oij/3xj  +  Fi  =  °»  eij  =  \  (3ui /3Xj  +  au j /a xi )  (18) 

c.j  =  [(1  +  v)  (Odo^)  -  v  5..{Ddokk}l/ER  (19) 

and  boundary  conditions, 

0 i j n j  =  Ti  on  ST  (20) 


A  V.  C-V  -V 


Observe  that  (19)  is  the  constitutive  equation  for  a  viscoelastic  material 
with  a  constant  Poisson's  ratio.  The  factor  Ep"^  is  introduced  to  cancel  the 
factor  Er  in  the  definition  of  {  }  in  (4). 

That  the  viscoelastic  solution  obtained  using  CP-III,  (1)  and  (2), 
satisfies  the  governing  equations  (18)  -  (21)  may  be  shown  by  substituting  it 
into  these  equations,  and  then  employing  (6)  -  (9)  and  the  fact  that  the 
elastic  solution  satisfies  (10)  -  (17)  and  that  3^/at  =  0  on  Sy;  the  last 
condition  excludes  rotation  of  the  normal  on  Sy  surfaces  which,  for  example, 
could  arise  from  melting  or  ablation.  It  should  be  noted  that  the  traction 

p 

T.,  (16),  in  the  elastic  problem  is  in  general  different  from  that  in  the 
viscoelastic  problem.  However,  the  condition  dSy/dt  <  0  stated  in  CP-III 

D 

implies  the  latter  is  known  at  all  past  times  and  thus  T.  may  be  explicitly 
calculated. 

When  considering  the  problem  of  bonding  of  two  different  materials,  a 
complication  arises  because  the  interfacial  stresses  acting  across  the  bonding 
surface  are  not  necessarily  equal  in  the  elastic  problem.  With  reference  to 
Fig.  2,  let  us  use  superscripts  (+)  and  (-)  to  denote  quantities  above  and 
below  the  surface  y=0,  respectively.  For  the  viscoelastic  problem  tT  +  TT  =  0 
along  this  surface,  whereas  from  (16)  the  elastic  traction  sum  is 

T*+  -  tJ"  =  {D+dT+  }  *  { D-dT 7  } 

=  { (D+-  D')dT+  )  (22) 

If 

Ej  D+  =  E‘  O'  (23) 


R+  tR- 
then  T.  +  T. 


=  0  or,  equivalently,  the  interfacial  elastic  stress  above  the 
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bonding  surface  equals  those  below  it.  If  D"  is  proportional  to  D+  we  can 
select  ER/ER  so  that  the  traction  sum  in  (22)  vanishes.  When  the  materials 
are  not  the  same  and  the  Poisson's  ratios  are  not  one-half,  normal  surface 
tractions  induce  material -dependent  shear  tractions  (unless  the  friction 
coefficient  vanishes)  and  thus  a  pure  mode  I  state  does  not  exist  at  the  tip 
[27]. 

In  order  to  avoid  mathematical  complexities  that  would  make  it  difficult 
to  understand  the  basic  physics  of  the  viscoelastic  bonding  problem,  we  shall 
assume  that  if  bonding  is  not  between  bodies  of  the  same  material  then,  at 
least,  their  Poisson's  ratios  are  one-half  and  their  creep  compliances  are 
proportional,  (23).  If  one  body  is  much  stiffer  than  the  other,  so  that  its 
deformation  can  be  neglected,  clearly  we  may  use  the  results  from  an  analysis 
based  on  these  conditions  even  if  the  relatively  rigid  body  does  not  satisfy 
them. 

4.  Results  from  crack-tip  analysis 

The  opening  displacement  v  above  the  bond  surface  in  the  neighborhood  of 
the  crack  tip,  Fig.  2,  may  be  obtained  from  CP- 1 II  and  the  elastic  solution  in 
[4-1 ), 

v  =  CR  fB  ojju1)  FU'/Odt'  (24) 

where 

FU7 0  5  \  [2U/t')1/2  -  in  [IJliL'Z&lj^l]  (25) 

*  1  -  U'/O  /z 

and 

CR  (1  -  v2)/Er  (26) 

The  surface  normal  stress  may  be  treated  as  a  specified  quantity,  and 

D 

therefore  its  counterpart  in  the  elastic  problem  is  obtained  from  (16), 


.  T/. 


r.  »  m 


II 


ft 


V. 


v* 

V 

v* 

>* 


which  yields, 

°b  =  {Ddob}  =  eR  JT  dT  {27) 

—  OD 

The  superscript^)  is  omitted  until  needed  later  to  distinguish  between 
quantities  above  and  below  the  bond  surface. 

The  stress  singularity  at  the  crack  tip  or  bond  edge  is  removed  by  using 
the  Barenblatt  result  for  elastic  media  141, 

KI  =  (f)1/2  Aj(Or1/2  <U  (28) 

o 

D 

where  Kj  is  the  stress  intensity  factor  for  an  elastic  material  with  Young's 
modulus  ER  and  Poisson’s  ratio  v.  This  stress  intensity  factor  is  found  from 

p 

an  elastic  analysis  for  which  ob  z  0,  and  is  considered  to  be  a  known  function 

p 

of  time  in  the  present  analysis.  Equation  (28)  implies  that  e  and/or  ob  must 

p 

change  with  Kj  if  the  stress  at  the  crack  tip  is  to  be  finite. 

5.  Prediction  of  quasi  steady-state  bonding  speed  when  ob  is  constant 

Equations  (24)  -  (28)  will  be  used  along  with  an  energy  criterion  for 

00  R 

bonding  to  predict  e  and  bonding  speed  ab  (=  -a)  as  functions  of  Kj  .  It  is 
assumed  that  ob  is  independent  of  time  and  is  spacially  uniform  over 
0  <  &  <  8.  Also,  ab  is  assumed  to  be  essentially  independent  of  time  during 
the  time  it  takes  for  the  tip  P  to  move  the  distance  s;  ab  may  vary 
appreciably  with  time  over  propagation  distances  which  are  large  compared  to 
e,  and  thus  the  terminology  quasi  steady-state  bonding  is  used.  In  the 
Appendix  a  formulation  involving  a  spacially  varying  ob  is  analyzed,  as 
mentioned  previously. 

Consider  now  the  process  of  bonding,  a  >  0,  and  identify  in  Fig.  2  a 
material  particle  (x  =  constant,  y  =  0)  which  is  currently  on  the  unloaded 
surface,  5>s.  Suppose  at  time  t=tj  the  left  end  of  the  bonding  zone. 
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5=6,  arrives.  This  particle  is  acted  on  by  the  constant  until  t=t£  when 
the  right  end  of  the  zone,  5=0,  arrives.  The  work  per  unit  surface  area  done 
on  the  particle  by  is  o^,  where  vfa  is  v  at  5=6,  assuming  that  v=0  at 
5=0.  The  lower  material  may  not  be  the  same  as  the  upper  material  and 
therefore  the  displacements  may  be  different.  This  difference  will  be 
identified  later;  but,  for  now,  we  shall  consider  only  the  displacement  above 
the  bond  plane. 

The  problem  at  hand  is  to  use  (24)-(28)  to  obtain  v^  and  thus  0^.  By 
adding  this  work  to  that  for  tne  lower  surface  and  equating  the  result  to  the 

so-called  bond  energy  denoted  by  2rfa  (per  unit  area)  an  equation  for 

§  R 

predicting  a&  is  obtained.  First  we  need  0^  ,  which  may  be  found  from 
(27)  by  noting  that  =  0  for  t<tj  and  is  constant  for  t^  <  t  <  t2',  thus, 
for  a  particle  at  fixed  x  and  with  tj  <  t  <  t£, 

%  -  er  %  <?9> 

The  relationship  between  t  and  5  follows  from  the  quasi  steady-state 
assumption, 

t  =  (t2  -  t)ab  and  8  =  (t?  -  t1)afc)  (30) 

from  which 

t  -  tj  =  (s  -  5)/ab  (31) 

Generalized  power  law  creep  compliance 

The  creep  compliance  of  many  materials  can  be  well -represented  over 
several  decades  of  time  by  the  so-called  generalized  power  law, 

D(t)  =  D0  +  (32) 

where  D0,  Dj  and  m  are  positive  constants  and  0  <  m  <  1.  We  shall  use  this 
form  to  obtain  exact  results  and  then,  guided  by  these  findings,  develop  an 


£ 


VJ 


16 


t 


§ 

i 


I 


s: 

** 


approximate  formula  for  more  general  forms  of  creep  compliance.  From  (29), 
(31),  and  (32), 

°b  =  ER  °blDo  +  (0/ab)mi  (33) 

where 

n  e  US  (34) 

Substitute  (33)  into  (28)  and  find 

kJ  =  (8/,)1/2  Er  ob  61/2  [D0  +  D1  yn  (e/ab)m]  (35) 

where 

i  (*/4)1/2  r(l*n)/r(1.5+m)  (36) 


I 


-\ 


and  :(•)  is  the  Gamma  function.  Observe  that  the  relationship  between  e  and 

R 

Kj  depends  on  bonding  speed,  whereas  the  analogous  relationship  for  crack 
growth  (which  relates  "failure  zone"  or  "fracture  process  zone"  length  a  to 
Kj)  does  not  depend  explicitly  on  speed  14-1].  Next,  use  (33)  in  (24)  to  find 
the  displacement  at  5=3,  and  thus  the  work  ;>bv  .  It  is  helpful  to  integrate- 
by-parts  the  term  involving  the  logarithm  and  obtain,  finally. 


:  Vb 


-?<!• 


2\  ^ 


eiD„  +  D,  c 


1mm 


(e/a.  )r 


(37) 


where 

cm  e  (2m+l)/(m*l)  (38) 

Previously  we  introduced  superscripts  (+)  and  (-)  to  denote  quantities 
above  and  below  the  bond  plane.  In  view  of  (23), 

=erd;  erd1  =er°;  <39> 

Recall  that  the  Poisson's  ratios  are  assumed  equal  and  constant.  The 
superscripts  can  be  omitted  in  (35)  as  only  these  products  appear.  However, 
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the  work  above  the  bond  plane,  which  will  now  be  denoted  by  W*,  may  be 
different  from  that  below,  W”,  since  ER  does  not  appear  as  a  factor.  Indeed, 
from  (37), 

/*£  •  °o  /»:  ■  i°\  <4o> 

The  criterion  for  bonding  may  be  written  as 

2rb  -  wj  .  U'  -  (1  +  C-  /0 l  )nj  (41) 

or,  equivalently. 


where 


s  v<>  ♦  °o  K  >  <43> 

For  identical  materials  (42)  reduces  to  =  W*  ,  whereas  if  the  lower 
material  is  rigid,  2rb  =  W*  .  We  shall  use  the  effective  bond  energy  in 
all  subsequent  work  so  that  the  notation  in  the  results  does  not  have  to  be 
changed  when  considering  the  bonding  of  materials  with  the  same  or  different 
stiffnesses;  the  effect  of  the  lower  material  will  appear  only  in  (43)  until 
we  apply  the  theory  in  Section  7. 

Equations  (35)  and  (42),  where  (37)  is  used  for  W^,  may  be  solved  to 

•  R 

obtain  e  and  ab  in  terms  of  Kj  ,  given  and  .  We  assume  these  latter  two 
quantitites  are  constants  and  obtain,  finally, 

B/Bm  =  Z(l+x)-2  cffl2  (44) 

where 

5  ’rb  «°b2<1-»2>°;  Sn2  <45> 
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is  a  constant  with  the  dimension  of  length.  Also, 


5  /5  .  Z(i«r2  C 

b  in  m 


where 


m  m  1  m  o  m' 


1/m 


(46) 

(47) 


is  a  constant  with  dimensions  of  velocity.  The  quantity  x  is  a  function  of  2, 
as  found  from  the  solution  to  a  quadratic  equation. 


x  =  0.5 (c  Z  +  (cmZ)1/2(cZ  +  4/c  -  4) 1/2 1-1 
m  '  m  m  m  ' 


In  turn, 


2  h  (kJ  /kJq)2 


(48) 

(49) 


where 


(O2  5  s )2/(‘-“2> 


(50) 


is  a  constant  with  dimensions  of  stress  intensity  factor  squared. 

For  the  range  1  <  Z  <  ®,  we  find  both  e  and  ab  are  monotone  decreasing 
functions  of  Z  and  that 


x  -  0  as  Z  -*  1  and  x  =  c  Z  for  Z  >>  1 

m 


(51) 


This  behavior  yields  from  (44), 


?  -1 
6/e  -  c  as  Z  -  1  and  6/ *  Z  for  Z  >>  1 
mm  m 


(52) 


and  from  (46), 


ab/am  -  «  as  Z  *  1  and  ab/am  =  Z"(1+1/m)  for  Z  »  1  (53) 


Figures  3  and  4  show  the  normalized  bonding-zone  length  (44)  and  bonding 
speed  (46)  as  functions  of  Z  on  logarithmic  scales  (log  e  log^).  The  range 
of  values  c  m  and  Z  for  which  the  behavior  is  essentially  that  for  the  high 
stress  inte  '  cy  factor  range  in  (52)  and  (53)  may  be  seen.  In  this  range  the 
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results  are  the  same  as  the  exact  solution  over  0  <  K j  <  ®  for  a  pure  power- 
law  material,  (32)  with  D0  =  0. 

In  order  to  obtain  explicit  analytical  results  for  Z  »  1,  em  and  am  may 
be  substituted  into  (44)  and  (46)  in  order  to  eliminate  DQ.  The  results  for  Z 
»  1  are 


R  ,2 


8  -  2.1 r*  Ej  /(l-v^tjcj  ) 


;(2+1/m)[(l_v2)C+  Ym/cmi1/m{Ej  /(l-v2)Kj)2(1+1/m)/8ob2cn 


(54) 

(55) 


These  forms,  i  .e. , 

B  =  k^K^/Ej)'2,  ab  =  k(K^/Ej)"q,  q  e  2(l+l/m)  (56) 

(where  k |  and  k  are  positive  constants)  are  like  those  for  the  crack  growth 
problem;  but  the  exponents  on  stress  intensity  factor  in  the  fracture  problem 

are  positive.  Here,  both  the  bonding-zone  length  and  speed  increase  with 

R  R  R  • 

decreasing  Kj.  Immediate  bonding  occurs  when  Kj  =  Kjq  since  ab  =  at  this 

2 

limit;  the  corresponding  value  for  e  is  6  c  ,  (52). 

R  R 

For  Kj  <  Kj  ,  (48)  yields  the  result  x<0.  A  physical  interpretation  of 

R  r  R  • 

this  low  Kj  behavior  may  be  made  by  noting  that  for  Kj  >  Kjq,  increases 

R 

with  increasing  effective  bond  energy  ,  for  Kj  fixed  (cf.  ( 45 ) - (50) ) .  When 
R  R 

Kj  <  Kjq  there  is  more  work  available  than  required  to  draw  the  surfaces 

R  R 

together,  neglecting  the  effects  of  material  inertia.  Hence,  for  Kj  <  a 
dynamic  analysis  would  be  needed  to  predict  finite,  high  speed  bonding. 


General  creep  compliance 

The  foregoing  results  based  on  the  generalized  power  law  compliance  (32) 
are  helpful  in  extending  the  theory  to  a  more  general  creep  compliance,  which 
we  write  in  the  form 
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D(t)  =  D0  +  aD( t )  (57) 

where  0Q,  as  before,  is  the  initial  compliance;  thus,  aD(0+)  =  0.  This 
compliance  and  (29)  -  (31)  yield 


°b  =  °bU)  =  ER  °b  fDo  +  ADt(8-0/*ab)}  (58) 

This  function  is  to  be  used  in  (28),  which  will  be  evaluated 
approximately  using  a  method  like  that  used  in  crack  growth  analyses  [4-11]. 
First,  rewrite  (28)  by  introducing  a  logarithmic  transformation  to  change  the 
integration  variable  to  l, 

L  e  1 og ( 1  -  r)  (59) 


where,  as  before,  log  e  logj0  and  r  =  Equation  (28)  becomes 


(-)1/2(in  10)s1/2 


°b  wldL 


where 


(i-ioLr1/2  ioL 


(60) 

(61) 


is  a  weight  function;  it  is  drawn  in  Fig.  5.  Because  this  function  is  very 
small  except  for  L  close  to  zero  (corresponding  to  e.=0),  and  aD  in  (58) 
decreases  with  decreasing  L  (or  increasing  0,  only  the  behavior  of  aD ( t )  for 
t  =  e/ab  needs  to  be  considered  (cf.  (58)).  Let  us  assume  that  aD  may  be 
approximated  by  a  power  law  aD  =  Djtm  over  at  least  a  one-decade  time 
range  (say),  i.e.  e/10ab  <  t  <  B/ab  (which  corresponds  to  -1  <  L  <  0)  and  that 
behavior  for  smaller  times  is  unimportant  because  w^  is  so  small  when 
L  <  -1,  (35)  may  be  used,  in  which 


m  = 


d  log  aD 
d  log  t 


(62) 


'*  /*  ;  ;• ,%  . 


t  Vs*  -  *  .  *  -«  *  -  \ 


where  t  =  s/ab. 


where 


A  more  compact  way  of  writing  (35)  is 


-  <!>1/2£Rv1/2|1(t'> 


*■  ■ 


(63) 

(64) 


It  is  found  that  t  is  practically  constant;  this  factor  is  a  monotone 

m 

increasing  function  of  m,  in  which  0.541  <  s  2/3  for  0  <  m  <  1. 

The  displacement  (24)  is  needed  to  obtain  the  work  where 

v&  ;  v  at  t=e.  Using  (59)  we  find 


where 


4  r  .  ...  r°  R 

vb  =  7  CR=(£n  10)  J  °b  w2  dL 
**  8 

w2  i  0.5  [2'_1/2  -  )10L 


(65) 

(66) 


and  -  =  1  - 1 0 L .  The  graph  of  the  weight  function  wp  in  Fig.  5  shows  that,  just 

D 

as  for  W| ,  the  behavior  of  ob  (and  thus  & D)  is  relatively  unimportant  for  L  < 
-1.  Just  as  (35)  has  been  replaced  by  (63),  we  may  replace  (37)  by 


wb 

=  *  (l-v2)ob26  D(t“) 

(67) 

where 

t" 

.'1/m  ft/a 
;  Ym  6/ab 

(68) 

and 

Ym 

r  cmYm 

(69) 

The  factor  is 

a  monotone  decreasing  function  of  m,  in  which 

1  i  Y' 1/m  <  1.47  for  1 
m 

>  m  >  0. 

The  bonding  criterion  (42)  may  now  be  written  as 


2rb  *  7  (‘-W8  °V> 


(70) 


'  v  •*- *  v  V  v  *c2x* itutm  >-V-  A'a.V.  Wivivilv...  v.  .W-V-  /■  A X.'/.  /■ 
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which  together  with  (63)  (recalling  that  E^D  =  E^D+)  provides  the  pair  of 
equations  for  finding  e  and  a^.  The  previous  results,  (44)  and  (46),  could  be 
used  as  first  approximations  in  an  iterative  solution  approach.  Here,  we 
shall  discuss  only  the  low  and  high  speed  limiting  conditions  for  bonding. 

Previously,  it  was  found  that  6  and  ab  are  monotone  decreasing  functions 

p  D  R  0  D 

of  Kj,  for  Kj  >  Kjq.  However,  ab  may  vanish  at  a  finite  value  of  Kj,  as  will 

be  shown  by  using  (70)  together  with  (54)  for  e.  Although  the  latter  equation 
is  exact  only  for  a  pure  power  law  compliance,  it  is  expected  to  be  a  good 
approximation  whenever  the  compliance  can  be  approximated  by  a  power  law  over 
one-decade  time  intervals.  (Observe  that  s  depends  on  m  through  cm,  (38),  but 
not  on  Di;  m  is  now  given  by  (62).)  Substitute  (54)  for  e  into  (70)  and 
obta  ,  i 

(<J)2  =  A  :■  0+(r)(Ej  )2/(l-v2)cm2  (71) 

R  4- 

where  t"  is  in  (68).  With  increasing  Kj,  D  (t")  must  increase  according  to 
(71).  However,  if  the  creep  compliance  has  a  finite  upper  limiting  value 
D+,  i.e.  D+  ;  D+U)  <  =  ,  then  the  largest  value  of  KR  for  which  a  solution 

CL  CD  ^  1 

a.  exists  is  K.  ,  where 
b  1® 

IK?.’*  5  S  o;(£j)2/(l-v2)  (72) 

When  t"  -  »,  then  D+  -  o+,  and  the  slope  m  vanishes;  thus  c„  -  1  as  t"  ♦  ®. 

A  further  study  of  (71),  in  which  it  is  assumed  that  m  in  (62)  does  not 

increase  with  t,  along  with  the  previous  results  (44)  and  (46),  leads  to  the 

R  R  R 

following  inequalities  and  behavior.  For  KJo<  Kj  <  K^,  the  length  e  and 

speed  ab  are  monotone  decreasing  functions  of  KjR;  also, 

e  -  b_c  ^  and  a.  -  ®  as  K?  -  KR  (731 

mm  b  i  lo 

a  -  8mc  ^  0 _/D  and  aK  -  0  as  K?  ♦  K?  (74) 

mm  o  ®  b  li® 
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Recall  that  D  /D  above  the  bond  surface  is  the  same  as  tnat  below.  From  (50) 
o  ® 

and  (72), 

-  «y°0)1/2  <«> 

6.  Comparison  of  apparent  fracture  and  bond  energies 

In  many  situations  involving  slow  loading  or  unloading  of  rubber,  the 
global  behavior  is  elastic,  with  this  behavior  defined  in  terms  of  elastic 
constants  which  are  essentially  equal  to  the  long-time  values  of  the 
associated  viscoelastic  properties.  On  the  other  hand,  at  this  same  time 
viscoelastic  effects  may  be  very  pronounced  in  the  neighborhood  of  a  moving 
crack  tip  or  contact,  edge  because  of  the  high  local  strain  rates. 

This  situation  is  easily  represented  in  the  theory  at  hand  by  taking  as 
the  aroitrary  reference  elastic  modulus  the  reciprical  of  the  actual  long-time 
compliance;  thus  D*  =  1/E^  and  0”  =  l/E”.  With  this  global  elastic  behavior, 
it  is  corranon  practice  to  define  an  apparent  fracture  energy  asociated  with 
crack  growth  by  means  of  an  elastic-like  equation  for  fracture;  in  terms  of 
the  notation  and  theory  at  hand,  in  which  adhesive  and  cohesive  fracture  are 
considered , 

■;f  M1-v2)D*  k|/4  (76) 

where  Kj  is  the  stress  intensity  factor  which  produces  a  particular  crack 
speed,  a.  Similarly  ,  guided  by  (72),  we  may  define  a  apparent  bond  energy 
using  the  same  expression, 

5  (I-”2)CT  <Ki>274  <”> 

R  • 

where  Kj  produces  a  bonding  rate  a^;  since  the  actual  global  behavior  is  also 

D 

that  of  the  reference  elastic  problem,  Kj  =  Kj. 


r, 

» 
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The  apparent  energies  in  (76)  and  (77)  can  be  related  to  the  (intrinsic) 
effective  fracture  energy  and  bond  energy  r^.  From  [4-IIJ,  after  allowing 
for  adhesive  fracture  through  an  equation  analogous  to  (43), 

r'  =  (l-v2)D+(a/3a)Kj/4  (78) 

and  by  eliminating  between  (76)  and  (78)  we  find 

raf=  rf  (79) 

The  compliance  ratio  is  the  same  above  and  below  the  bond  surface,  as  assumed 
previously.  Similarly,  from  (71)  and  (77), 

r;b  =  r'  D(f)/D,cn2  (80) 

The  ratio  D^/D ( t )  increases  monotonical ly  (from  unity)  as  time  t  decreases 
from  infinity  to  zero,  and  for  rubber  D  /D„  may  exceed  one  hundred  [28], 
Hence,  the  ratios  r^/r^.  and  rt/rab  may  themselves  exceed  one  hundred. 
Moreover,  if  the  effective  energies  are  equal,  =  r^,  the  ratio  of  apparent 
energies  is  bounded  according  to 

1  5  —  5  ('mW2  <8I> 

rab 

In  experiments  on  the  time-dependent  adhesion  of  rubber  to  glass,  Roberts 
and  Thomas  (111  found  that  the  ratio  r^f/r^b  1ndeed  ver>'  Targe  over  a  wide 
range  of  bonding  and  debonding  speeds.  This  ratio  was  reported  to  exceed  one 
thousand  at  the  shortest  testing  time,  which  was  approximately  ten  seconds. 

Finally,  it  is  of  interest  to  observe  for  this  case  of  global  elastic 
behavior  that  (72)  reduces  to 

■  4rb/0l(1-v?>  (82) 


3 


a  v,  A  w  /.  V,  *■.  •*.  •*,  A  *. 


css 
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and  that  ab  >  0  if  Kj  <  K^.  For  the  fracture  problem  a  >  0  if  Kj  >  Kj^, 


where 

(kJj2  2  4r;/0*(l-.2)  (83) 

Consequently,  if  the  bond  and  fracture  energies  are  equal,  (82)  or  (83)  define 
a  value  of  the  stress  intensity  factor  below  which  there  is  bonding  and  above 
which  there  is  crack  growth. 

7.  Examples 

Cohesive  crack  healing 

Consider  the  Griffith  type  of  problem  in  which  there  is  an  isolated 
circular  crack  with  radius  a,  or  a  through-the-thickness  crack  of  length  2a. 
A  remote,  uniform  tensile  stress  o(t)  normal  to  the  crack  plane  is  given.  In 
the  reference  elastic  body  the  stress  intensity  factor  is  |29|, 

kj  =  c  a1/2oR  (84) 

1/2  1/2 

where  c  is  a  constant;  c  =  2/n  for  a  circular  crack  and  c  =  n  for  a 

R 

through-crack.  Equation  (16)  or  the  inverse  of  (1)  provides  o  , 

oR  =  Er  J  D(t-x)^  dx  (85) 

—  oo 

which  is  assumed  to  be  non-negative  for  (84)  to  be  valid.  Let  us  use  the 
power  law  in  (56)  for  ab;  then  with  (84)  for  KR  and  with  ab  =  -a,  we  find 

a/aQ  -  (1-It)1/P  (86) 


where  a0  is  the  size  at  t=0.  Also, 


p  =  q/2  +  1  =  (l+2m)/m 


! 
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\  -=  *<2pa”p  Jt(coR/ER)’qdt  (88) 

o 

The  healing  time  th  is  defined  as  the  time  at  which  a=0;  namely,  It(th)  = 

An  explicit  solution  for  th  may  be  easily  found  if  o  is  a  constant  which  is 

D 

applied  long  before  t=0,  so  that  (85)  yields  o  =  EdD  o,  where  D  e  D(°>)  is 
the  equilibrium  or  long-time  creep  compliance.  Equation  (88)  reduces  to 


and  thus 


lt  =  k2?ao’P  (cDa.°^"q  1 

*h  =  *0P(cMq/k2P 


(89) 

(90) 


Observe  that,  as  expected,  the  higher  the  tensile  stress  o,  the  longer  the 
healing  time.  Also,  notice  that  the  arbitrary  constant  ER  does  not  appear  in 

p 

the  prediction  of  a  and  th  because  o  /ER  is  independent  of  ER  (cf.  (85)). 

It  should  be  recalled  that  whether  or  not  healing  or  crack  growth  occurs 

R  R 

depends  on  the  value  of  Kj  relative  to  in  (72).  For  the  cohesive  healing 

problem,  =  r b/2  and  0*  =  D^.  Using  also  (84)  and  (85)  the  criterion  for 

bond  growth  is  found  as 

ca1/2  J  D(t-0  dx  <  l2rbDm/(l-v2)]1/2  (91) 

-  cc 

For  a  pure  power  law  compliance  (which  includes  that  for  a  Newtonian  body, 
m=l)  =  «,  and  (91)  is  always  satisfied  if  (o,t)  <  ®.  But  crack  growth 
analysis  shows  that  the  criterion  for  growth  also  is  always  satisfied  (4- 
II J.  This  quasi-steady  state  analysis  is  thus  not  sufficient  to  determine  the 
sign  of  a  for  a  power  law  material. 

Growth  of  contact  area  between  initially  curved  surfaces 
Figure  6  shows  the  contact  region  for  two  bodies  which  are  pressed 
together  by  the  load  F  and  are  also  drawn  together  under  bonding  stress  o^. 


r.- 


The  length  scale  e  over  which  acts  is  assumed  to  be  very  small  compared  to 
the  contact  size.  Just  as  in  the  fracture  problem  (i.e.  contact  area 
reduction)  the  opening  displacement  is  parabolic  near  the  contact  edge  except 
very  close  to  the  edge  where  the  cusp  shape  of  Fig.  2  exists,  shown  in  Fig.  6 
by  the  blackened  edge  area.  In  order  to  predict  afa(=  ac)  the  stress  intensity 

D 

factor  for  the  reference  elastic  problem  Kj  is  needed.  Here,  we  shall 
consider  two  cases:  (i)  contact  between  initially  spherical  surfaces  with 
radii  R+  and  R"  and  (ii)  the  same  type  of  problem,  except  the  surfaces  are 
cylindrical  so  that  the  deformation  is  two-dimensional.  An  approach  to 

D 

finding  Kj  will  be  used  which  is  similar  to  that  in  (12)  for  contact  between  a 
plane  and  a  spherical  surface. 

First,  in  the  absence  ofa.  Hertz's  solutions  for  contact  size  may  be  used 
for  the  spherical  surfaces, 

ac3=  3F*/8BE^  (92) 

and  the  cylindrical  surfaces, 

ar2=  2F*/„BE^  (93) 

P 

where  FH  in  (92)  is  the  total  force,  while  in  (93)  it  is  force  per  unit  length 
along  the  cylinder  axis.  Also, 

B  =  ( 1/R+  +  l/R')/2,  ~  s  (-“  +  —)  (1-v2)  (94) 

F 1  F  F~ 

tR  lR  lR 

d 

Now,  let  us  reduce  the  compressive  load  by  an  amount  aF  ;  assume  the 
surfaces  are  bonded  so  that  ac  is  fixed  at  the  Hertz  value.  This  change 
produces  an  interface  tensile  stress  which  is  infinite  at  the  contact  edges 
(without  the  Barenblatt  modification),  because  these  edges  are  really  crack 
tips.  The  stress  intensity  factor  may  be  obtained  from  the  solution  for  the 


i< 
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normal  stress  between  a  bonded  flat-ended  rigid  punch  under  tensile  load  aFk 
and  flat  elastic  halfspace  {251,  using  the  assumption  that  ac«R+  and  ac<<R“; 
the  Hertz  solutions,  (92)  and  (93),  are  based  on  this  assumption,  which  we 


retain.  That  this  is  the  same  normal  stress  as  produced  by  aFk  in  the  problem 


in  Fig.  6  may  be  argued  by  first  noting  that  it  is  independent  of  elastic 


constants  and  the  initial  surface  radii,  and  produces  a  rigid  translation  of 


the  contact  surface.  The  same  normal  stress  acts  on  both  upper  and  lower 


surfaces,  and  since  the  contact  surface  shape  does  not  change  the  upper  and 


lower  bodies  remain  in  contact.  Because  the  interface  conditions  are 


satisfied  using  this  normal  stress  and  the  bodies  are,  in-effect,  halfspaces 
(i.e.  ac<<R+  and  ac«R“),  it  is  the  one  and  only  solution  for  interface 


stress.  (However,  the  conditions  for  zero  interfacial  shear  stress  must  be 


met,  as  discussed  at  the  end  of  Section  3.) 


For  a  circular  punch  with  radius  ar,  the  contact  stress  is, 


5  a  /J  _.2n  1/2 
2i,ac(ac  -  J  } 
aFk 


2*ac(ac+  r)1/2(ac-  r)1/2 


where  r  is  the  radial  distance  from  the  contact  center  to  a  point  on  the 


interface.  From  the  singularity  at  r=ac  and  the  definition  of  stress 


intensity  factor. 


K j  =  11m  [2n(ac-  r)l1/2oJ 


we  obtain 


K*  -  AFR/2„1/2ac3/2 


Similarly,  for  the  plane  strain  problem, 


■as 


m 


29 


kJ  =  AFR/(,ac)1/2  (98) 

Now,  aFR  =  FR  -  FR,  where  FR  is  the  total  load.  With  this  expression 

D 

and  Fh  from  (92)  and  (93),  the  stress  intensity  factor  in  (97)  for  spherical 


surfaces  becomes. 


K*  -  48E;a  3/2/3.‘'2  -  ,  -3«f«/2.1/2 


For  cylindrical  surfaces. 


K*  ■  „1/2BE‘a3/2  n  -  a;1/2FRA1/2 


(100) 


p 

The  force  FR  in  the  elastic  problem  is  related  to  the  actual  compressive  force 

R  R 

through  F  =  (DdFj,  just  as  the  stress  in  (85).  Introduction  of  Kj  in  terms 

of  a£  makes  (99)  and  (100)  differential  equations  for  ac. 

As  simple  examples,  the  contact  si2e  will  be  found  analytically  from 

these  differential  equations  using  the  power  law  bonding  speed  equation  (56) 

for  two  cases:  FK  =  0  and  Kj  =  0.  In  the  first  case,  (99)  yields  the  radius 


ac  '  A  V 


(101) 


where  t^  is  the  bonding  time,  and 


s  5  2/(3q+2) 


A  =  (k(3*1/2D'/4B)q/s)S 


D'  r  (1  +  DJ/D|)(1-V2) 


(102) 


(103) 


(104) 


The  last  equality  comes  from  (39).  In  deriving  (101)  we  assumed  interfacial 
contact  starts  at  t=0  and  neglected  the  error  in  the  Barenblatt  method  at 
short  times  (when  ac  <  e).  The  size  ac  for  contact  between  cylinders  has  the 
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same  time  dependence  as  in  (101)  because  (100)  depends  on  ac  in  the  same  way 
as  (99)  when  =  0.  From  (56)  and  (102), 

s  =  m/(3+4m)  (105) 

If  m=l  (viscous  material)  then  s=l/7,  which  is  the  largest  physically 
acceptable  exponent. 

The  time  dependence  of  ac  in  (101)  is  in  general  considerably  different 
from  that  predicted  from  the  Hertz  contact  problem.  Taking  for  example  a 
constant  load  F  applied  at  t=0,  the  load  in  the  elastic  problem  is 
FR  =  F^D^tbmF.  From  (99)  and  (100),  respectively. 


ac  =  (3D'D*  F/88) 1/3  tbm/3 

(106) 

ac  =  (20'D*  F/nB)1/2  tbm/2 

(107) 

(Anand's  model  (131  predicts  this  time-dependence  for  a  power  law  material.) 

□ 

Now,  let  us  retain  Kj  but  assume  it  is  small  enough  that  ac  is  close  to 
the  Hertz  solutions  (106)  and  (107).  For  the  power  law  of  (56),  this  occurs 
when  ab  is  sufficiently  large.  Let  aH  be  the  Hertz  solution,  which  is  for 

p 

Kj  =  0  in  (99)  and  (100).  Then  ac  =  +  aa,  where  aa/aH«  1  by  assumption. 

Using  this  approximation  in  (99)  we  find 

aa/aH  =  i.1/20‘  (k/aH)1/q/4BaH3/2  (108) 

p 

where  aH  is  the  solution  to  (99)  for  Kj  i  0.  For  example,  using  (106), 
aa/aH  -  tb2m  /3(m+l) ,  anc)  thus  the  8rror  in  the  Hertz  solution  (106) 
diminishes  with  time.  A  similar  result,  with  the  same  exponents  as  in  (108), 
is  found  for  the  contact  between  cylinders;  using  (107)  ,  we  find  again  that 
the  error  diminishes  with  time  in  that  aa/aH  -  . 
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Effect  of  bonding  time  on  joint  strength 

As  discussed  in  Section  1,  the  time-dependent  strength  of  bonds  formed 
across  contacting  surfaces  has  received  considerable  attention.  The  examples 
considered  in  the  present  section  deal  only  with  contact  area  growth  and  not 
the  processes  following  contact,  such  as  polymer-polymer  interdiffusion  (which 
may  be  very  significant  if  the  temperature  is  sufficiently  high).  Here,  we 
observe  only  that  (99)  and  (100)  may  be  used  to  predict  strength,  if  the 
residual  viscoelastic  effect  of  the  bonding  problem  is  negligible  so  that  the 
correspondence  principle  for  crack  growth  may  be  used.  In  this  case,  the 
elastic  and  viscoelastic  stresses,  loads, and  stress  intensity  factors  are  the 


same.  Suppose,  for  example,  that  following  the  bonding  process  the 

temperature  is  reduced  enough  that  the  critical  stress  intensity  factor  Kj^  is 
a  constant  and  that  the  first  term  in  (99)  is  negligible  at  fracture. 
Then  the  strength  is  F  -  a^2  KjC.  Using  (101), 

F  -  while  from  (106)  F  -  t^2;  for  viscous  media  m=l,  and  the 

21  5 

strengths  are,  respectively,  F  -  t^'  and  F  ~  tfa*  . 

If,  following  bonding,  we  use  the  idealization  from  [20)  in  which  any 

nonuniformity  in  stresses  in  the  bond  is  ignored,  use  elastic  fracture 

mechanics,  and  assume  the  fracture  energy  r^.  is  proportional  to  bonded  area 

1/21/2 

Afa,  then  KjC  ~  ~  Ab  '  ;  with  a  given  crack  at  the  interface  the  joint 

1/2  s/2 

strength  is  proportional  to  Afa  .  From  (101),  KIC  ~  tb  '  for  spheres  and 
KIC  ~  V  for  cylinders;  if  m=l,  these  exponents  are  1/14  and  1/7, 
respectively.  On  the  other  hand  from  (106)  and  (107),  -  tbm^  and 

KIC  -  respectively. 

For  polymer-polymer  bonding  of  noncrossl inked  systems  somewhat  above  the 

1/4 

glass  transition  temperature,  the  behavior  -  tb  is  usually  found 
experimentally  [19,20]  (in  a  range  for  which  m  =  1  is  probably  a  good 
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assumption);  it  is  also  the  behavior  predicted  from  interdiffusion  theory  for 
fully  contacted  or  wetted  surfaces.  None  of  the  exponents  discussed  above 
yield  the  1/4  exponent  (for  the  physically  acceptable  range  of 
viscoelastic  exponents  0  <  m  <  1)  except  for  the  last  one  of  m/4,  if  m=l,  for 
contacting  cylindrical  surfaces.  Anand  [17)  reported  that  when  nearly  flat 
surfaces  are  bonded,  the  initially  irregular  contacting  areas  become  in  time 
contacting  cylinders,  and  that  the  conversion  to  cylindrical  contact  surfaces 
occurs  early  in  the  bonding  process.  For  this  case,  one  could  not  use  an 
experimentally  determined  exponent  of  1/4  to  conclude  that  the  increase  in 
joint  strength  is  due  to  interdiffusion  as  opposed  to  contact  growth.  Of 
course,  if  the  externally  applied  contact  pressure  is  high  enough  during 
bonding  so  that  essentially  full  contact  is  quickly  established,  contact  area 
growth  will  not  be  an  important  source  of  time-dependence  of  the  joint 
strength. 
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8.  Conclusions 

A  mechanics-based  theory  of  bonding  between  the  same  or  different  linear, 
isotropic  viscoelastic  media  has  been  developed  and  used  in  some  examples.  We 
assumed  quasi  steady-state  bonding  and  the  bonding  stress  c>b  and  energy  to 
be  independent  of  bonding  speed  ab  in  order  to  obtain  explicit,  closed-form 
solutions.  However,  if  these  material-related  quantities  vary  with  speed 
(which  could  be  due  to  surface  roughness  on  a  scale  which  is  very  small 
compared  to  the  bonding-zone  length  b)  most  of  the  equations  do  not  have  to  be 
changed.  For  example,  with  speed  dependence  of  these  parameters,  (55) 
becomes  an  implicit  equation  for  ab;  with  power  law  dependence,  the  resulting 

9  D 

solutions  for  e  and  ab  obey  power  laws  in  Kj,  as  In  (56),  but  the  exponents 
are  in  general  different. 
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Regardless  of  how  or  if  and  depend  on  bonding  speed,  we  may 
conclude  that  all  dependence  of  ab  and  6  on  external  loading  and  geometry  of 
the  continua  is  through  the  instantaneous  stress  intensity  factor  Kj.  This 
simplicity  exists  as  long  as  8  is  small  compared  to  local  geometric  scales 
such  as  initial  radii  of  surface  asperities  or  crack  length;  as  mentioned 
above,  very  small-scale  roughness  compared  to  6  is  also  acceptable. 

p 

In  contrast  to  the  crack  growth  problem,  Kj  depends  on  external  loading 

•  R 

history,  and  afa  and  e  decrease  with  increasing  Kj.  The  external  loading  must 
be  such  that  afa  is  nonnegative  in  order  for  the  correspondence  principle  used 
in  the  analysis  to  be  valid.  The  correspondence  principles  used  for  bond 
growth  and  crack  growth  are  different.  If  crack  healing  occurs  after  crack 
growth,  the  bonding  theory  will  not  be  valid  unless  effects  of  the  prior 
history  of  crack  growth  have  essentially  faded  out.  The  complexity  of  the 
analysis  appears  to  be  much  greater  when  the  velocity  of  the  boundary  of 
bonded  surfaces  changes  sign  and  viscoelastic  effects  from  before  and  after 
the  change  are  combined;  however,  methods  exist  for  analyzing  this  type  of 
problem  [ 30,31  ] . 

Linear  theory  may  not  be  valid  under  high  stresses,  such  as  when  bodies 
with  initially  curved  surfaces  are  pressed  together  under  high  compressive 
loads.  In  such  cases,  a  J  integral  theory  may  be  applicable.  It  would  be 
similar  to  that  for  crack  growth,  but  based  on  correspondence  principle  III 
(instead  of  II,  which  is  for  crack  growth)  given  in  [26]. 

Although  the  development  here  is  for  isotropic  media,  through  a  change  in 
the  creep  compliance  D(t)  the  theory  in  Sections  5  and  6  can  be  used  for 
bonding  of  surfaces  of  orthotropic  media  if  the  surfaces  and  bond  edge  are 
parallel  to  principal  material  planes.  This  extension  is  based  on  the 
findings  in  [5]  for  crack  growth,  In  which  the  equations  for  the  mechanical 


*  JV 

•Cv 


w  m  w  vmmrnmrnmmmnmmmmmmmmmim 


34 

state  near  and  at  the  interface  are  the  same  as  for  isotropic  media;  in  this 
case  (1-v  )  D(t)  is  to  be  replaced  by  another  function  of  time  consisting  of  a 
combination  of  creep  compliances  of  the  orthotropic  material.  (The  same 
generalization  can  be  made  in  applications  involving  effectively  unbounded 
media,  such  as  those  in  Section  7.)  However,  just  as  for  isotropic  materials, 
property-dependent  interfacial  shear  stresses  develop  unless  the  bonding  is 
between  the  same  materials  or,  if  different,  between  incompressible  materials 
with  proportional  compliances;  one  material  may  be  rigid.  The  present  theory 
has  been  developed  under  the  condition  that  there  is  no  shear  stress  at  the 
interface  in  order  to  bring  out  the  essential  features  of  the  mechanics  of 
bonding  or  healing  of  viscoelastic  materials  with  a  minimum  of  mathematical 
complexity. 

Finally,  it  should  be  noted  that  Poisson's  ratio  appears  in  elastic 
p  2 

\  solutions  in  this  paper  only  in  the  form  (1-v  )/modulus,  as  in  (26)  and 

(94).  For  problems  in  which  this  is  the  case,  one  can  show  by  means  of 

» 

'  Graham's  extended  correspondence  principle  [61  that  the  theory  remains  valid 

|  when  Poisson's  ratio  is  time-dependent,  as  long  as  it  is  the  same  for  the 

materials  above  and  below  the  bond  surface.  In  this  situation  one  may  use  a 

,  p 

plane-strain  creep  compliance  C(t)  in  place  of  (1-v  )  D(t)  throughout  the 

paper;  this  replacement  is  analogous  to  the  generalization  made  for 

orthotropic  media,  as  discussed  above. 
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APPENDIX 


The  stress  in  the  bonding  zone  was  assumed  constant  in  Section  5. 


Here  we  shall  use  the  form 


°b  =  °of(v) 


(A.  1 ) 


where  v  =  (v+  -  v*)/z0,  in  which  oQ  and  zQ  are  constants  with  dimensions  of 
stress  and  length,  respectively.  The  separation  between  the  surfaces  is  v+- 


v".  The  bond  energy  is 


2rh  =  z  o  J  f  dv 

D  0  0  J 


Greenwood  and  Johnson  [12],  in  a  study  of  adhesive  crack  growth  used 


f  =  (1  +  v)-3 


(A. 2) 


(A. 3) 


based  on  experimental  and  theoretical  considerations  of  the  force  of 
attraction  between  smooth  surfaces.  Figure  7  shows  their  predicted  distri¬ 
butions  of  f  and  v  for  elastic  materials,  where  Ep  is  defined  in  (94).  Here 
we  will  not  use  (A. 3);  the  only  assumption  about  f  we  make  is  that  it  vanishes 
rapidly  enough  with  increasing  v  to  ensure  validity  of  the  Barenblatt  method 
(i.e.,  a  small-scale  effective  bonding-zone  length)  and  convergence  of  the 
relevent  integrals.  In  addition,  we  shall  use  the  creep  compliance  (32),  but 
with  D0  =  0,  and  assume  as  before  that  the  compliance  of  the  material  below 
the  bonding  surface  D"  is  proportional  that  above  this  surface,  D+;  hence. 


e;°;  =er°; 


(A.4) 


The  surface  stress  in  the  elastic  problem,  (27),  for  quasi  steady-state 


bonding,  (31),  becomes 


sr*' 


(A. 5) 
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with  df/dt1  =  df/dv  •  dv/dt'.  From  (24)  and  (A. 4)  with  s  =  “, 

2  CT> 

V  =  (v+  -  V)/z0  =  (1  +  D“/D|)  fo^U'^U'/Odt'  (A. 6) 

ER  Zo  0 

Substitute  (A. 5)  into  (A. 6)  and  eliminate  t  and  t'  in  favor  of  the 
dimensionless  variable 

*  =  i'U)  =  (Di  +  O'  )oc  tm+1/abniz0(m+l)  (A. 7) 


Also,  define  if'  i  y(s’)  and  v' '  i  hence. 


v  = 


0  * 


/  <l> 1  '  \ p  i  im  df 

(  ’  ‘  11  d*1 


F[U7*)p]d*'  1  d  ii>  ’ 


(A. 8) 


Here,  p  i  l/(m+l).  Inasmuch  as  f  -  f(v),  this  is  an  integral  equation  for 
v  =  v(i>).  The  primary  result  for  our  purposes  is  that  v  is  a  function  of  only 
the  dimensionless  variable  *.  This  observation  together  with  (A. 5)  permits  us 
to  reduce  (28)  to 


kJ/E 


+ 

R 


k2zo 


/2^y  5/C*1) 


(A. 9) 


where  t<2  =  l<2(m)  is  dimensionless  function  of  only  m.  By  solving  for  ab  the 
same  exponent  q  as  in  (56)  is  obtained.  In  order  to  determine  the  value  of  k 
in  the  power  law  (56)  it  would  of  course  be  necessary  to  find  k£,  which 
requires  the  solution  of  (A. 8).  It  should  be  added  that  Greenwood  and 
Johnson  [12)  obtained  a  ~  in  the  crack  growth  problem  where  q  is  the  same 
as  in  (56). 

We  may  use  arguments  similar  to  those  in  Section  5  to  show  that  (A. 9)  is 
valid  for  a  compliance  which  does  not  obey  a  pure  power  law,  as  long  as 
d  log  0/d  log  t  is  a  slowly  varying  function  of  log  t.  However,  without 
further  analysis  it  is  not  possible  to  indicate  just  how  slow  the  variation 
must  be  to  achieve  an  acceptable  degree  of  accuracy. 
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Figure  2.  Normal  stress  and  displacement  along  the  bond  surface  in  the 
neighborhood  of  the  bond  edge. 
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